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ABSTRACT 

We have used the E p ^ - Ei so correlation of GRBs to measure the cosmological param- 
eter J7m- By adopting a maximum likelihood approach which allows us to correctly 
quantify the extrinsic (i.e. non-Poissonian) scatter of the correlation, we constrain 
(for a flat universe) SIm to 0.04-0.40 (68% confidence level), with a best fit value of 
Om ~ 0.15, and exclude flu = 1 at > 99.9% confidence level. If we release the as- 
sumption of a flat universe, we still find evidence for a low value of Hm (0.04-0.50 
at 68% confidence level) and a weak dependence of the dispersion of the E p ^ - Ei so 
correlation on f^A (with an upper limit of J7a ~ 1-15 at 90% confidence level). Our 
approach makes no assumptions on the E p ^ - E[ so correlation and it does not use 
other calibrators to set the "zero' point of the relation, therefore our treatment of the 
data is not affected by circularity and the results are independent of those derived 
via type la SNe (or other cosmological probes). Unlike other multi-parameters corre- 
lations, our analysis grounds on only two parameters, then including a larger number 
(a factor ~ 3) of GRBs and being less affected by systematics. Simulations based on 
realistic extrapolations of ongoing (and future) GRB experiments (e.g., Swift, Konus- 
Wind, GLAST) show that: i) the uncertainties on cosmological parameters can be 
significantly decreased; ii) future data will allow us to get clues on the "dark energy" 
evolution. 

Key words: gamma-rays: observations - gamma-rays: bursts - cosmology: cosmo- 
logical parameters. 



1 INTRODUCTION 

Gamma-ray Bursts (GRBs) are the brightest cosmological 
sources in the universe, with isotropic luminosities up to 
10 54 erg cm -2 s _1 and a redshift distribution extending at 
least up to z ~6.3 (e.g., Tagliaferri et al. 2005). Thus, these 
sources may be interesting for cosmological studies, if one 
can use them to provide measurements of the cosmolog- 
ical parameters independently of other methods, like the 
CMB (e.g., De Bernardis et al. 2000; Spergel et al. 2003; 
Dunkley et al. 2008 ; Komatsu et al. 2008), type la SN 
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(e.g., Perlmutter et al. 1998; Perlmutter et al. 1999; Riess 
et al. 1998; Riess et al. 2004; Astier et al. 2006), Baryon 
Acoustic Oscillations (BAO) (e.g., Eisenstein et al. 2005; 
Tegmark et al. 2006; Percival et al. 2007), galaxy clusters 
(e.g., Rosati, Borgani & Norman 2002; Voit 2005). However, 
GRBs are not standard candles, given that their luminosi- 
ties span several orders of magnitude under the assumption 
of both isotropic and collimated emission ( e.g., Ghirlanda 
et al. 2006). In the recent years, several attempts to "stan- 
dardize" GRB have been made, mainly on the basis of cor- 
relations involving: a) a GRB intensity indicator, like the 
isotropic-equivalent radiated energy (Ei so ) or the isotropic- 
equivalent peak luminosity (L p ,i so ); b) the photon energy 
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at which the time averaged vF„ spectrum peaks ("peak en- 
ergy", E Pt {); c) other observables, like the break time of the 
aftergl ow light curve, tb, or the "high signal time scale" 
Tp. 45 |Ghirlanda, Ghisellini fc Lazzatil |2004| ; iFirmani et al.l 
l2006l ). Finally. ISchaefeij l|2007h used different correlations 
(spectrum-energy, time lag - luminosity, peak luminosity - 
variability) to construct a GRB Hubble diagram and con- 
strain the cosmological parameters. 

These analyses have provided useful constraints on Qm 
and Qa consistent with those derived from type la SNe, 
see, e.g., iGhirlanda et al. I (|2006l ) for a review. However, the 
use of these correlations for cosmology is controversial. For 
example, because of the lack of low redshift GRBs they 
cannot be directly calibrated. On the other hand the cal- 
ibration of a spectrum-energy c orrelation using S Ne-Ia as 
calibrators (|Kodama et al. 2008; Liang et al.ll2008l ) is an in- 
teresting attempt to use GRBs as cosmological probes, but 
it may look suspicious that this method gives very similar 
results to those obtained via SNe. Thus, particular care and 
sophisticated statistical methodologies have to be adopted 
in order to avoid circularity problems when constructing a 
GRB Hubble diagram. In addition, recent analyses based 
on updated samples of GRBs showed that the dispersion 
of three-parameter c orrelations could be significantly larger 
than thought before (ICampana et al ]|2007l ; IGhirlanda et al. I 
120071 : iRossi et al.ll2008l ) 

In this article we investigate the possibility of con- 
straining the cosmological parameters from the E Pi i - 
-Ei so correlation. This correlation was initially discovered 
on a small sample of_BeppoSAX GRBs with known red- 
shift (|Amati et al.l l2002h and confirmed afterwards by 
Swift observations l|Amatil 120061 ). Although it was the first 
"spectrum-energy" correlation discovered for GRBs (and 
the most firmly established one, to date) it was never used 
for cosmology purposes, because of its significant "extrinsic" 
scatter (i.e., a scatter in excess to the "intrinsic" Poissonian 
fluctuations of the data). However, the conspicuous increas- 
ing of GRB discoveries combined with the fact that E Pi i - 
Ei so correlation needs only two parameters that are directly 
inferred from observations (this fact minimizes the effects of 
systematics and increases the number of GRBs that can be 
used, by a factor ~ 3) makes the use of this correlation an 
interesting tool for cosmology. 

Our study was motivated by our finding that, in the 
assumption of a flat universe, the \ 2 obtained by fitting the 
E P] i - Ei so correlation with a power-law varies with the value 
of Qm assumed to compute Ei so , has its minimum value 
for J7m~0.25, very close the value obtained with SNe-Ia ap- 
proach. In this paper we show indeed that a small fraction 
of the extrinsic scatter is due to to the choice of cosmologi- 
cal parameters, thus allowing us to constrain Qm and, to a 
smaller extent, Qa- 
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Eiso SAMPLE 



Previous analyses of the E Pt i - Ei so plane of GRBs shows 
that different classes of GRBs exhibit different behaviours: 
while normal long GRBs and X-Ray Flashes (XRF, i.e. 
particularly soft bursts) follow the E Pi i - Ei so correlation, 
short GRBs and the peculiar very close and sub-energetic 
GRB 980425 do not (Amati et al. 2007). Therefore, in our 



analysis we considered only long GRB/XRF. The sample of 
long GRB/XRF used in this paper (updated to April, 11 
2008) is reported in Table 1 and includes 70 GRBs. We re- 
port, for each GRB, the redshift, the fluence and the cosmo- 
logical rest-frame spectral peak energy, £" p ,i=iJp,obsX (1+z). 
The redshift distribution covers a broad range of z, from 
0.033 to 6.3, thus extending far beyond that of Type la 
SNe (z <~1.7). For Swift GRBs, the E p ,i value derived 
from BAT spectral analysis alone were cons ervatively taken 
from t he results reported by the BAT team l|Sakamoto et alj 
200830)- Other BAT E p>i values reported in the literature 
we re not considered , becaus e either they were not confirmed 
by ISakamoto et "all (2008a, b) refined analysis (e.g., Cabr- 
era et al. 2007) or they are based on speculative methods 
(|Butler et al.ll2007h . 

In Table 1 we also report the values of Ei so computed 
from published spectral parameters and flue nces by follow- 
ing the standard method described, e.g., by lAmatl (|2006l ) 
and assuming Ho = 70 km s _1 Mpc -1 , Qm = 0.3 and 
SI a = 0.7 (references are reported in Table 1). An analy- 
sis of Figure 1 (left panel) shows that, for this "standard" 
cosmology, all GRBs in our sample follow the E Pt i - Ei so cor- 
relation; in particular, Swift GRBs are well consistent with 
the E Pi i - Eiso correlation based on GRBs discovered by 
other instruments characterised by different trigger thresh- 
olds and sensitivities (as a function of energy). This fact 
points out that the E p ,i - Ei so correlation is not affected by 
significant selection effects (see also Ghirlanda et al. 2008). 

If we fit the data of Figure 1 with a simple power- 
law, we find an index m = 0.57 ± 0.01 and a normalization 
K = 94±2, consistent with results of previous analyses (e.g., 
Amati 2006). However, despite the very high significance of 
the correlation (Spearman's p = 0.872 for 70 events), the fit 
with a power-law provides an highly unacceptable x^ vam e 
(408/68). This is clear evidence of the existence of a signifi- 
cant extrinsic scatter, which is superimposed to the Poisso- 
nian one and implies the existence of "hidden" parameters, 
connected with GRB phenomenology, which contribute to 
define the location of a GRB in the E p ^ - Ei BO plane. It 
should be noted that systematics affecting both E Pt i and -Eiso 
may contribute to the extrinsic scatter, even if there is evi- 
dence that they are not a dominant component (|Landi et al.l 
l2006h . 
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Eiso CORRELATION 



The main goal of this paper is to investigate if the extrinsic 
dispersion of the E Pi i - Ei so correlation is sensitive at varying 
the values of the cosmological parameters Qm and Qa- We 
emphasize that this method does not suffer from circularity, 
since we do not assume an E Pt i - Ei so relation based on a 
particular choice of the cosmological parameters or calibrate 
it by using other cosmological probes. 

As a first step, under the assumption of a flat universe, 
we studied the trend of the x 2 values obtained by fitting 
the E Pt i - Eiso correlation with a simple power-law, for 
different choices of Qm- Here and in the the following, we 
assumed Ho = 70 km s -1 Mpc -1 . Figure 2 (left) shows 
that the extrinsic scatter of the correlation indeed decreases 
with cosmology, and minimizes for Qm ~0.25. This result is 



Cosmology with the E p ^ - E [so correlation in GRB 

Table 1. Values of redshift, "bolometric" fluence and cosmological rest-frame spectral peak energy, -B p ,i=-E Pl obs X (l+ z )i 
of long GRBs and XRFs with firm estimates of both z and E p ^ s (70 events) as of April, 11 2008. These are the values 
that we used to estimate cosmological parameters. The table includes also the values of E ISO computed by assuming 
Ho =70 km s — 1 Mpc - , Sly = 0.3 and Ha = 0.7 . The uncertainties are at lcr significance. The "Instruments" column 
reports the name of the experiment(s), or of the satellite(s), that provided the estimates of spectral parameters and 
fluence (GRO = CGRO/BATSE, SAX = BeppoSAX, HET = HETE-2, KW = Konus-Wind, SWI = Swift). The last 
column reports the references for the spectral parameters used to derive the fluence and E p ^. Redshift values were taken 
from the GRB table by J. Greiner and references therein (http://www.mpe.mpg.de~jcg/grbgen.html ). 
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Notes. ( a ) Bolometric fluence computed in the [1/(1 + z) - 10000/(1 + z)\ keV energy range. 

( b ) Computed by assuming Ho =70 km s -1 Mpc -1 , !!m = 0.3 and Q\ = 0.7 . 

( c ) References for the spec tral paramete r s and for the values and uncertainties of E Pt i and Ej sn : (11 [ Amatil ||2006|) 
and references ther ein; (21 lUlan ov et a hl J2005T) and refined analysis of - BejPjPoSAXdata; (3) ISakamoto e t al.l J2 005t) ; (4) 
IStratta et all [|2006ll : (5llCenko et al.l ll2006h~7ollKrimm et al.N2006h : (7l lSakamoto et alJ J2u08ah : ( 81 iGolenetskii et al. I 
ll2006ar); (9) ICampana et al.l ll2006h ; (10) IGolenetskii et al. I (120061 ) ; (111 ISchaeferl ||2007|) and r e ferenc es therein; (12) 
lAmati et al 1 &200H) and references t herein; (13) IGolenetskii et al I J2006J|: (141 iMundell et al.l ll2007l) and references 
ther ein; ( 151 iGhirlanda. et al. I ||2007|1 and references therein; ( 161 iPerlev et al.l ll2008l); (171 IGolenetskii et a"I~l ll2007ah; 
(181 IGolenetskii et al. I j2007bll:*7l9l IGolenetskii et al. I ll2007dl: ( 201 IGolenetskii et al. I l|2007dh : (211 IGolenetskii et al. I 
^OOSall ; (221 IGolenetskii et al I ll2008bh ; (231 IGolenetskii et al I l l2008dl . 




Figure 1. Left: location in the E p ^ - E- lao plane of the 70 GRBs and XRFs with firm estimates of redshift and B Pi0 b s included in our 
sample (Ei so computed following Amati (2006) and assuming a cosmology with Ho = 70 km s — 1 Mpc -1 , Om = 0.3 and £7a = 0.7). Red 
dots are Swift GRBs. Black dots are GRBs discovered by other satellites. The best— fit power— law is the continuous line (±2<r GX t region). 
Right: distribution of the normalised scatter (see, e.g., Rossi et al. 2008 for definition and method) of the _E p ; - Ei so correlation for 
Hq = 70 km s — 1 Mpc -1 , Qm = 0-3 and = 0.7); the normalised Gaussian is superimposed to the data. 



in qu alitative agreement with the one obtained, e.g ., with 
SNe JPerlmutter et all Il998l , Il999l ; iRiess et all Il998l . 12004 
lAstier et al.ll2006l ). 

In order to better characterize the dependence of the 
_Ep,i - Ei BO extrinsic scatter on cosmology, we have used 
the maxi mum likelihood m ethod (hereafter MLM) as dis- 
cussed bv lD'Agostinil (|2005l ) (see also Guidorzi et al. 2006, 
Amati 2006 and Rossi et al. 2008 for other applications of 
this methodology). This method assumes that, if the cor- 



related data (x;,y;) can be described by a linear function 
Y — mX + q with the addition of an extrinsic scatter <T cxt 
among the free parameters, the optimal values of the param- 
eters (m, q, and cr cxt ) can be obtained by minimizing the - 
log(likelihood) function, in which the uncertainties on both 
(x;,yi) are taken into account. The general log(likelihood) 
function is given by: 

logp(m,q,<7 x ,a y \{xi,y t ,a x ,i,ay : i}) = 
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Figure 3. Value of the extrinsic scatter <r ex t (left) and power— law index m (right) as a function of SJm, obtained by fitting the correlation 
with the MLM method which accounts for extrinsic variance, for a flat universe. In the left panel, we also show, as a dashed line, the 
result obtained by using a different log(likelihood) (see text). 



1 N 1 

2 Yl [ l0g ( 2 n(a* + m 2 a 2 +a 2 + 

(yi-mxj-q) 2 1 
a 2 + m 2 a 2 + a 2 yi + m 2 a 2 .^ J 

where, in our case, log(-E Pl i), o x — and a y = cr cxt . We re- 
mark that setting cr x = does not affect the general validity 
of this formula. 

As an example, the fit of the correlation with this 
method by using the Ei so values reported in Table 1 provides 
the following parameters value: m — 0.54 ±0.03, K — 98 ±7 
and o"oxt = 0.17 ± 0.02 (68% c.L; consistent with the results 
reported by Amati 2006 for a smaller sample of 41 GRBs). In 
Figure 1 (right) we show the distribution of the normalised 
scatter of the data (see, e.g., Rossi et al. 2008 for details 
on the computation of this quantity) with superimposed the 



normalised Gaussian. It is apparent that the scatter of the 
E Pt i - Eiso correlation is indeed dominated by the extrinsic 
(i.e. non-Poissonian) variance. 

For the goals of this paper, we repeated the above anal- 
ysis by varying !5m, always under the assumption of a flat 
universe. Figure 2 (right) and Figure 3 (left) show that the 
values of the -log(likelihood) and of the extrinsic variance 
cr cx t are indeed sensitive to Qm, both showing a clear min- 
imum around Qm ~0.15. Also the slope of the correlation 
is sensitive to the adopted cosmology, as shown in Figure 3 
(right). By using the probability density function, the MLM 
also allows us to constrain S)m, that results to be in the 
range 0.04-0.40 at 68% and in the range 0.02-0.68 at 90% 
c.l (Table 2). An Qm value of 1 can be exluded at ~99.9% 
c.L 

If we release the flat universe hypothesis and let Qm 
and Qa vary independently (Figure 4, left), we still find ev- 
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Table 2. 68% and 90% c.l. ranges of Qm (both by assuming a 
flat universe or by letting Cm and Qa to vary independently) as 
determined with different methods applied to the present sample 
of 70 GRBs. 



Method 


c.l. 




(flat) 




n A 


scatter 


68% 


0.04 


- 0.40 


0.04 - 0.50 


<1.05 




90% 


0.02 


- 0.68 


0.01 - 0.75 


<1.15 


scatter 


68% 


0.04 


- 0.40 


0.05 - 0.41 


<1.05 


(self calib.) 


90% 


0.02 


- 0.67 


0.01 - 0.73 


<1.13 


scatter 


68% 


0.03 


- 0.28 


0.03 - 0.33 


<1.10 


(m = 0.5) 


90% 


0.01 


- 0.49 


0.01 - 0.53 


<1.17 



idence for an universe with a low value of Qm (0.04-0.50 
at 68% c.l.). Only an upper limit of ~1.05 can be set to 
Qa (see Table 2). This fact is not surprising, given that 
the redshift distribution of GRBs is expected to produce 
very vertically elongated contours in the Q m - Qa plane 
|Ghirlanda et al. IbOOfj ; iKodama et~aI1l200St ). The comple- 
mentarity of GRB with other cosmological probes can be 
seen in Figure 4, where we plot, in addition to our res ults, 
the contours derived from SN la bv lAstier et al. | <|2006h . As 
can be seen, already with the present sample of GRB, sig- 
nificantly improved constraints on Qm can be obtained by 
the combination of the two probes. 

We have also investigated the estimate of the free para- 
maters of the the inverse relation X — mY + q plus the 
extrinsic scatter parameter a ex t ~ cr x - The results are fully 
consistent with previous ones and satisfy the expectation 
that \frnm' = r xy , where r xy is the Pearson's weighted 
linear-correlation coefficient (e.g., Bevington 1969, Bendat 
& Piersol 2000). For a flat universe with Qm= 0.3, we find 
m = 0.54, m = 1.57, Vmm' = 0.92 and r xy = 0.92. 

Finally, we t ested our r esults using the likelihood func- 
tion suggested by iReichartl (|200ll ). The results obtained are 
fully consistent with the results above (Figure 2, left). This 
shows that the constraints that we have derived on the cos- 
mological parameters do not depend significantly on the 
adopted statistical methodology and confirms the general 
soundness of our approach. 



4 FUTURE PERSPECTIVES 

4.1 Enrichment and redshift extension of the 
sample 

In order to estimate the accuracy on cosmological parame- 
ters achievable by ongoing and future experiments, we have 
carried out a number of simulations based on an enriched 
sample of GRBs. By using the slope, normalization and ex- 
trinsic scatter of the E Pt i - Ei so correlation obtained with 
the real sample of 70 GRBs and assuming, as an example, 
ilM=0.27, fiA=0.73, we generated 150 random GRBs with 
an accuracy on E Pi i and Ei so of 20%. For 90% of the sim- 
ulated GRBs, the redshift was randomly chosen from the 
GRB redshift distribution of our sample. For the remain- 
ing 10% of simulated GRBs we assumed a value of redshift 
larger than 6. The number of GRBs, and the accuracy on 
E Pt i and Ei so , of the simulated sample are those that will re- 
alistically be provided in the next years by joint observations 
of Swift and GLAST/GBM (plus also Konus-Wind, Suzaku, 



RHESSI, AGILE, etc.). Also the extension of the redshift 
distribution up to at least z = 9 is consistent both with the- 
oretical predictions based on star formation rate evolution 
and the expected improvement in the sensitivity of hard X- 
ray and optical/near-infrared telescopes, that will be used 
to detect and follow-up GRBs and measure GRB redshift 
(e.g., Salvaterra et al. 2007 Salvaterra & Chincarini 2007). 

In Figure 4 (right) we show the confidence level contours 
for finiand Qa obtained by applying the method described 
in Section 3 to the sample made of the 70 "real" GRBs plus 
the 150 "simulated" ones. In Table 3, we report the con- 
straints on cosmological parameters both in the assumption 
of a flat universe and by varying simultaneously Qm and 
Qa- A simple comparison with Table 2 shows that the sim- 
ulated sample decreases the uncertainty range of Qm from 
10 to only a factor 2. In addition, as can be seen in Figure 4 
(right), significantly improved constraints on both Qm and 
Qa will be obtained by the combination of the contours ob- 
tained from GRBs with those from other probes, like, e.g., 
SN la. 

Finally we note that the E Pi i - E-, correlation, because 
it needs the measure of a third observable (i.e. the break time 
in the optical afterglow light curve) requires the estimates 
of z and E p ,i for a number of GRBs larger by more than 
a factor of 3 than required by the E Pi i - Ej so correlation. 
For in stance, the simulations performed bv lGhirlanda et al. I 
(2006) with 150 GRBs correspond to a "real" sample of at 
least ~500 GRBs having both z and E Pi j measured. 

4.2 Calibrating the E Pi \ — Ei so correlation 

One of the main problems with the use of GRB correlations 
for cosmology is the lack of low-redshift GRBs (i.e. lying at 
z <<1) allowing for a cosmology-independent calibration 
similarly to type la SNe. This problem can be partly over- 
come by fitting the correlation in a sub-sample of GRBs 
lying at similar redshifts. This technique was proposed for 
future calibration of the slope of the E Pt i - E-, correlation 
(e.g., Ghirlanda et al. 2006), but can be effectively used for 
the E P: i - Ei so correlation, due to the much larger number 
of events that can be included in the E Pt i - Ei so sample 
with respect to the E Pt i - E y one. By studying the redshift 
distribution of the GRBs in our sample, we find that there 
are 18 events lying at redshift between ~0.8 and 1.2. The 
analysis of this sub-sample with the MLM method provides, 
indeed, evidence of a very low variation of m as a function 
of cosmological parameters. For instance, by assuming a flat 
universe, we find that the m best-fit values are in the range 
m=0.493-0.496, to be compared to m=0.477-0.560 obtained 
on the whole sample of 70 GRB. The value of m holds sta- 
ble within this interval even by varying both Qm and Qa 
in the [0,1] range. The dependence of m on the dark en- 
ergy equation of state (parametrised as discussed in next 
section) is also within this interval. However, when taking 
into account the uncertainty resulting from the fits, the con- 
straint on the calibrated slope of the E Pt i - Ei so correlation 
becomes larger: m= 0.44-0.55, for both a flat and general 
universe. Thus, the applicability of the method is currently 
hampered by the large uncertainty affecting m, which re- 
sults in a marginal improvements of the c.l. of cosmological 
parameters (Figure 5, left; Table 2). On the other hand, if 
we apply this calibration method to the simulated sample of 
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Figure 4. Contour confidence levels of S!m and Q^, obtained by fitting the correlation with the MLM method which accounts for 
extrinsic variance, with the present sample of 70 GRBs (left) and with the improved sample of 70 + 150 GRBs expected from next 
experiments (right; see text). In both panels we also show as blue dashed contours the constraints derived from SN la by Astier et al. 
(2006). 




1.5 




Figure 5. Contour confidence levels of S!m and Ha obtained by fitting the correlation with MLM method which accounts for extrinsic 
variance with m calibrated based on GRBs lying at 0.8 < z < 1.2 (black contours) and m fixed at 0.5 (red contours); see text. Left: 
present sample of 70 GRB; right: improved sample of present 70 GRBs + 150 GRBs expected from next experiments. 



70 + 150 GRB, the improvement in the estimates of cosmo- 
logical parameters is significant (see Table 3 and compare 
Figure 5, right, with Figure 4, right). 

A second possibility to calibrate the slope of the cor- 
relation is to derive it on robust physical basis. Several 
analyses reported above or elsewhere (e.g., Ghirlanda et al. 
2008) point out a slope of m=0.5. This value is also pre- 
dicted by several theoretical models, involving in various 
forms synchrotron, inverse Compton, t hermal and Comp- 
tonise d thermal emissions; s ee, e.g., IZhang fc Meszaroa 
l|2002l) ; lThompson et all i|2007T ). Figure 5 (red contours) and 
Table 3 illustrate quite clearly that the correlation fitted 
with the MLM method after assuming m = 0.5 improves 
significantly the accuracy degree of the cosmological param- 
eter measurements. 



4.3 Investigating the evolution of the dark energy 

The GRB redshift distribution extends up to z >~ 6, which 
is well above that of SNe-Ia (z ~ 1.7). Therefore, at least 
in principle, GRBs are powerful tools to study the evolution 
of dark energy with the redshift. In the following, we adopt 
the parame trization of the dark energy equation of state 
propo sed bv lChevallier fc Polarskil l|200ll ); lLinder fc Uthererl 
|2005l ). i.e. P = w{z)p, where: 

w(z) = wo + ^f- z (2) 

As shown in Figure 6 (left), with our sample of 70 GRBs 
we find that the -log(likelihood) computed with the MLM 
by assuming a flat universe with Qm = 0.27 and by set- 
ting w a — is indeed sensitive to wq, with a minimum 
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Figure 6. Left: Value of the -log(likelihood) as a function of the parameter of the dark energy equation of state too ( see text) in the 
assumption of a flat universe with !Jm = 0.27 and Qa = 0.73 obtained with the present sample of 70 GRB. Right: contour confidence 
levels of £7m and too m the assumption of a flat universe obtained with the improved sample of present 70 GRBs + 150 GRBs expected 
from next experiments (continuous line: w a = 0; dashed lines: to a =4). In both figures, the red lines indicate the results obtained by fixing 
m at 0.5 . 



Table 3. 68% and 90% c.l. ranges of Qm (both by assuming a 
flat universe or by letting Qm and Qa to vary independently) as 
determined with different methods applied to a sample composed 
by the present sample of 70 GRBs plus 150 simulated GRBs with 
known z and E p j expected to be provided in the near future by 
present and planned GRB experiments (see text). 



Method c.l. (2m (flat) Qm 



scatter 


68% 


0.21 


- 0.50 


0.11 


- 0.57 


0.35 - 1.12 




90% 


0.16 


- 0.66 


0.08 


- 0.73 


<1.15 


scatter 


68% 


0.22 


- 0.48 


0.12 


- 0.49 


0.47 - 1.10 


(self calib.) 


90% 


0.17 


- 0.59 


0.10 


- 0.69 


0.10 - 1.11 


scatter 


68% 


0.13 


- 0.27 


0.18 


- 0.28 


0.62 - 1.10 


(m = 0.5) 


90% 


0.11 


- 0.35 


0.06 


- 0.35 


0.50 - 1.12 



around —1 (corresponding to a cosmological constant). Un- 
fortunately, within the limits of present GRB sample no sig- 
nificant constraints on wo can be provided. Even after as- 
suming m = 0.5, Wo can be constrained only to be between 
—3.7 and —0.7 at 68% c.l. In Figure 6 (right) we show the 
Qm - wo confidence level contours obtained with the sample 
of 70 real + 150 simulated GRBs by assuming a flat universe 
with w a — (continuous lines) and w a — 4 (dashed line). 
As can be seen, such a sample would provide clear evidence 
of Wo<— 1 and some hints on w a - 



5 CONCLUSIONS 

We have used an updated sample of 70 GRBs aimed at deriv- 
ing the cosmological parameters Hm and Qa from the E Pl i - 
Ei so correlation of GRBs. With respect to previous attempts 
mainly based on the E Pt i - E-y and E Pt i - £"i so -To.45 correla- 
tions, the use of E Pt i - Ei so correlation has the advantages of 
requiring only two parameters, both directly inferred from 
observations. This fact allows for the use of a richer sample 
of GRBs (e.g., a factor 3-4 larger than used in the E Pi i - 



E-y correlation), and a reduction of systematics, with respect 
to three-parameters spectrum-energy correlations. 

Our method consists of studying the dependence of the 
dispersion of the E p ^ - Ei so correlation on cosmological pa- 
rameters by adopting a maximum likelihood method, which 
allowed us to parameterise and quantify correctly the ex- 
trinsic (i.e. non-Poissonian) scatter component. From our 
analysis, a number of significant results do emerge: 

a) a simple power-law fit shows a clear trend of the 
dispersion, with a minimum at Qm = 0.2-0.4 (for a flat 
universe). By refining the fit with the MLM, we found that 
both the extrinsic variance and the slope of the correlation 
show a significant dependence on cosmology. From the study 
of the log(likelihood) we derived Q M = 0.04-0.40 at 68% c.l 
and 0.02-0.68% at 90% confidence levels; 

b) after assuming "a priori" a slope m = 0.5, consis- 
tently with both observations and the predictions of several 
models for GRB prompt emissions, we derive Qm = 0.01- 
0.49 at 90%; 

c) if we release the assumption of a flat universe, we still 
find evidence for a value of Q M < 1 (0.04-0.50 at 68% c.l.) 
and a weak dependence of the dispersion of the E Pt i - Ei so 
correlation on Qa with an upper limit of 1.15 (90% c.l.); 

d) our study does not make assumptions on the E Pt i - 
Ei ao correlation or make use of independent calibrators to 
set the "zero point" of the relation, therefore our approach 
does not suffer from circularity and provides independent 
evidence for the existence of a gravitationally repulsive en- 
ergy component ("dark energy") which accounts for a large 
fraction of the energy density of the universe; 

e) in a flat universe, the E Pz i - Ei so correlation is also 
(weakly) sensitive to the wo parameter of the equation of 
state of dark energy; 

f) we have simulated the impact of ongoing/planned 
GRB experiments (e.g., Swift + GLAST) on the future Qm 
and Qa measurements and shown that the uncertainties can 
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be decreased by almost an order of magnitude with respect 
to those obtained with the current GRB sample. 
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